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Abstract 

The problem of synthesizing a controller for plants subject to arbitrary, finite en¬ 
ergy disturbances and white noise disturbances via Linear Matrix Inequalities (LMIs) is 
presented. This is achieved by considering white noise disturbances as belonging to a 
constrained set in /. In the case of where only white noise disturbances are present, 
the procedure reduces to standard 34 2 synthesis. When arbitrary, finite energy dis¬ 
turbances are also present, the procedure may be used to synthesize general mixed 
performance objective controllers, and for certain cases, Robust 34 2 controllers. 


1 Introduction 

In the standard robust control paradigm, the signal space which characterizes perfor¬ 
mance is equivalent to that which captures a system’s uncertainty. For example, 34 x 
tools are used when dealing with bounded energy (or power) gain uncertainty (see [12], 
[10], [18]), while when working with disturbances, the uncertainty is assumed to be 
of finite amplitude gain (see [9]). While it is often the case that the particular charac¬ 
terization of the uncertainty is not critical to the design process, the signal space used 
to characterize the performance often is. In particular, one of the common complaints 
among control design engineers which use 344 methods is that the resulting designs 
tend to be sluggish and overly conservative. As an alternative, 3-4 2 designs are often em¬ 
ployed, but they lack the robustness properties of 34 M designs (see [5]) which can readily 
be extended to encompass a system’s uncertainty. The attractive feature of 34 2 designs 
is their gain interpretation; they minimize the power output when the disturbances are 
assumed to be white noise or impulses. This is in contrast to 344 designs, which mini¬ 
mize the energy to energy (or power to power) gain; in many applications, modeling the 
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disturbances as arbitrary signals is a poor modeling choice, and thus designs lead 
to low performance controllers. 

A desirable control design strategy would be one which has the input-output gain 
interpretation of the 3~C 2 norm, but can readily accommodate 3P X bounds on the uncer¬ 
tainty. 

In [15], a framework is developed whereby white noise signals are captured in a de¬ 
terministic setting. The main motivation behind this approach was the reconciling of 
the worst case setting, natural when considering robustness issues, with the stochastic 
setting. This framework proved very natural when addressing the so-called Robust 3~C 2 
Analysis problem, which was solved in [13] and [14]. 

This approach will be used in this paper to tackle the problem of Robust J-f 2 Synthe¬ 
sis for a restricted class of problems; in particular, rank one synthesis problems with 
time varying uncertainty. This will be achieved by solving an auxiliary problem, that of 
synthesizing a controller for plants subject to arbitrary, fini te energy disturbances and 
white noise disturbances. 

Other work in this area includes [2], [6], [19] and [17]. Also worth noting is research 
on the so called mixed J-C 2 / problem, where only nominal stability and no m i n al 3-f 2 
performance are considered (see [1], [8]). 

The paper is organized as follows: we begin with some mathematical preliminaries, 
followed by a review of the notions introduced in [15] with regards to capturing noise 
signals as elements of a set. The problem is then posed and solved, followed by a dis¬ 
cussion on computational issues. The problem of robust disturbance rejection is then 
addressed, and it is shown how these types of problems can be re-formulated in the gen¬ 
eral problem setup previously solved. We conclude with an example which illustrates 
the tools developed and their numerical properties. 

2 Preliminaries 

Most of the notation in this paper is standard. We restrict ourselves to discrete time sys¬ 
tems. The space of square summable sequences is denoted l 2 ; when the spatial structure 
is relevant, it is referred to as l v 2 . The 2-norm of a signal d in l 2 is denoted \\d\\. The 
unit ball of l 2 signals is denoted BI 2 , and consists of all l 2 signals whose norm is less 
than or equal to one. The discrete time, unit delay operator is denoted A. The truncation 
operator Pt is defined as 


(P r x)(t):= 


x(t) for t < T 
0 for t > T 


( 1 ) 


Causal, finite dimensional, linear time invariant systems will be denoted FDLTI. A causal 
linear map G over l 2 is bounded if the restriction of G to l 2 is a bounded operator, with 
induced l 2 norm denoted ||G||. The transfer function representation of FDLTI system G 
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is denoted G( A). The linear fractional transformation (LFT) between two systems G and 
K is denoted G * K, and is defined as: 

G * K\=Gn + G 12 K(I — G22K ) 1 G 2 \ (2) 


where 


G = 


G 11 G\ 2 
G 2 \ g 22 


when the inverse of (I - G22K) is well defined. BA is the set of linear, but otherwise 
arbitrary, operators whose induced 2-norm is less than or equal to 1. For two subsets of 
l 2 , Si and S 2 , the approximation error of Si with S 2 is defined as 


D(Si,S 2 )\= sup inf ||5i-5 2 || (3) 

5i e5i ^eS 2 

D(Si,S 2 ) is a measure of how well S 2 captures the elements of Si. For a constant matrix 
A, its maximal element is denoted 


||A|| m := max | A itj \ 
uj 


(4) 


3 Deterministic Noise Sets 


We begin by reviewing the notions introduced in [15] to capture white noise in sets. Given 
a signal nei™, its autocorrelation function is defined as 


00 

Rn(T):= X n(t)n r (t + T) 

t = — 00 


( 5 ) 


Note that there is no time averaging in the definition above, as would be used for power 
signals, since we are dealing with finite energy signals. Given positive integer N and 
positive number y, we define the following set of autocorrelation functions: 




R(-t) = R T (r) 

rntn — 
-tv N,y ~ ' 

R(t) : 7L — 

IH?(0)-/|| M <y 



11 JR(r) || M ^ Y, 1 < t < N 


and corresponding signal set 


= {n e lf\R n (T) e n% y } (7) 

In particular, when y = 0, we have 


VV N W Ni o 


( 8 ) 


It is shown in [15] how the worst case gain from IT™ to l 2 approaches the 'hfi norm 
in the limit as N goes to infinity, the rate of convergence being exponential. 
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4 Problem Formulation 



The problem formulation is as follows; given positive integer N and FDLTI system Go 


Go 


I 

II 

Cr n 

G12 

•l 

21 

rn 

Cj 2 i 

G22 


(9) 


find internally stabilizing (see [21]) FDLTI system K such that 

sup ||(G 0 * K)do\\ 2 < 1 (10) 

ZEBl 2 ,neW\f 

The performance objective is absorbed into plant Go; in general, the performance 
objective will be decreased until a controller no longer exists which satisfies (10). We 
will later show how the solution to the above problem may be used to solve a variety 
of robustness problems. The above problem formulation is interesting in its own right, 
however; it may be the case that some of the disturbance signals may be arbitrary, in 
this case l, while the others are white noise, such as n. For example, this may be the 
case when tracking a reference signal l (which may be weighted to restrict tracking over 
a certain frequency range), in the presence of sensor noise or other random disturbance 
n. 


5 Solution 

We will provide a solution to the above problem in a series of steps. The first consists 
of recalling the solution of controller synthesis when the disturbances are subject to 
implicit constraints, presented in [3], The next step consists of parametrizing sets TV,” 1 
in image form. The final step is to combine these image representations with the solution 
provided in the first step to solve (10). 
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5.1 Synthesis with implicit constraints 

The following is a review of the main results in [3]. Consider the following subset of l 2 




d el 2 


\\Hid\\ 2 < 1 

I Ml 1 2 = I Ml 1 2 


1 < i < c 
1 <j<c 


( 11 ) 


where the Hi, Lj and /, are constant matrices. Consider the following constrained fea¬ 
sibility problem: let system G and set J-f be given. Find a stabilizing controller K such 
that 


sup ||(C? * K)d\\ 2 < 1 (12) 

It is shown in [3] that such a controller K exists if and only if there exist sy mm etric 
matrices S, T, W, Z (of a given spatial structure) such that 


V R 1 


S 0 
0 I 


R 


S 0 

0 z 


V T < 0 


U T R 


T 0 
0 W 


R 1 


T 0 
0 I 


U < 0 


T I 
I S 


> 0 


z / 
I w 


> 0 


(13) 


C T ZC < 1 

where U, V, R and C are constant matrices which depend on the state space represen¬ 
tation for G and the constraint set . A state space representation for K may then be 
constructed from R, S and T (the details may be found in [3] and [11]). The above is a 
convex feasibility problem, and may be solved using numerical packages such as LMILab 
[7]. 
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5.2 Image Representation for IT™ 


The solution provided in the previous section cannot be utilized to directly solve (10); in 
order to directly specify IT™, Hi, Lj and Jj need to be systems, not constant matrices 
(see [2]). What we would like to do is construct an alternate representation for X\P which 
is consistent with the solution provided in the previous section. 

Let N and m be given. Define 


14(A) := 

V k ( A) := 

U := 
V := 


1 + A k 
V8N 
1 - A fc 
V8N 

[Ui Ui ... U N U n ] 
diag [U,U, U ] 


1 < k <N 


(14) 


where V e 2 , ie., V consists of m copies of U along the “diagonal”. Then it 
may be verified that 


VV* 



ie., V and U are co-inner, and \ \v\\ 

II 11oo 

Define n e l 2 ™^' as follows: 



i 

y2iv• 


n := (ni, ..., n m ) 
fli .= (fl-i ,i, ..., 71 


(15) 


(16) 


and the following set of constraints: 


Ci: 

<■ 7luk > 

7li,k) < 1, 

(ni,k,n itk > < 1 

1 < i < m, 

1 < k < N 

C 2 : 

(7li,k, 

7lj,k) ~ (n, 

i,fe> 7ij,k) = 0 

1 < i < j < m, 

1 <k<N 

C 3 : 

(7li,k, 

fij,k) - <n, 

i,k, 7lj,k) = 0 

1 < t < j < m. 

1 < k < N 

C 4 : 

(TLi'k, 

rij,k) + (il, 

i,k, 7lj,k) ~ 0 

1 < i < j < m, 

k = 1 


The constraint set 3\f is then defined as: 

-2V:={n e C2, C3, C4 are satisfied} (17) 

The image set AT™ y may now be defined: 

1 W™ y \={n E Z^ln = Vn,n E N, Until > 1 - y for 1 < i < m} (18) 

For y = 1, we define 

^n-=^Z (19) 

which corresponds to no explicit norm constraint on nt. The following theorem ou tlin es 
how Xv and Xv!y are related, and is crucial to the synthesis results which follow: 
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Theorem 1 


j /J w m 
W N, 0 


TV ; m 


AT 


2 . AT™ c 1 v^ y for y > 0 . 

3. Di'WJf'y, lY™) is upper semi-continuous as a function ofyaty-O. 
Before proving Theorem 1, we will need the following two Lemmas: 


Lemma 1 

S:={n e I™! 2NV*n e 2sf, llndl 2 = 1 for 1 < i < m} = iv ™ 0 (20) 


Proof of Lemma 1: 

It is clear by setting n = 2NV*n that S c AV™ 0 ; we will thus show that c S. Let 
n G T^ 0I with corresponding n. For each component, by constraints Ci, Hndl 2 < 2iV, 
which implies Hndl 2 < 1. Thus \\nt\\ 2 = 1, ||n<| | 2 = 2iV. 


Since t) is co-inner, 3U X e R.3-f (2N l)x2N such that V2N 


U 

A 

f/x 


is unitary. Thus for 


each i, Ui can be uniquely decomposed as 


n< = £/*V( + U*Wi (21) 

where v* e h, wt g 1 2N_1 - Furthermore, llndl 2 = IIII 2 + ||(/*Wj|| 2 - From this we 
may conclude that Vj = 2Nnu and that ||Vi|| = \\U*Vi\\ 2 = Undl 2 = 2N; thus Wj = 0, and 
nt = 2NU*rii. This gives n = 2NV*n, as required. 

□ 


Lemma 2 Given R e 21™ y , where y < —^ +1)2 , there exists a ^z'gfnai 1 x e ZJ such fhat 
J? x (t) = K(t) forr G [—AT,N]. 


The proof of Lemma 2 may be found in the Appendix. We are now in a position to 
prove Theorem 1: 

Proof of Theorem 1: 

We begin by showing how n is constrained when 2 NV*n g thf. Let n = 2NV*n. Thus 

n uk = 2 NU*m (22) 

ni,k = 2NUjfni for 1 < i < m, 1 < k < N 

For the above, it can be verified that constraints Cj to C 4 are equivalent to : 


Ci: 

lint 

II 2 + (nu\ k n.i) 

< 1 

1 

< 

i 

< 

m, 


1 

< 

k 

< 

N 


llw-i 

II 2 - (nu \ k rii) 

< 1 

1 

< 

i 

< 

m, 


1 

< 

k 

< 

N 

C 2 : 

(nt, 

, A k rij) + (nu A 

l 

s 

II 

0 

1 

< 

i 

< 

j * 

m, 

1 

< 

k 

< 

N 

C 3 : 

(nu 

A k nj) - (nu A 

0 

II 

s 

1 

< 

i 

< 

j < 

m, 

1 

< 

k 

< 

N 

C 4 : 

(nu 

nj) = 0 


1 

< 

i 

< 

j ^ 

m 
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Proof of 1: Let n e 'W™. The above constraints are then trivially satisfied, proving 
n e S; thus by Lemma 1, n e 'W™ 0 - Now let n e TT™ 0 , or equivalently, n e 5. It is then 
straightforward to show that the above constraints imply that n e lY™, as required. 

Proof of 2: This follows from 1. and c ir™ yo for 0 < y < yo- 


Proof of 3: Let e > 0 be given. It will be shown that there exists a yo > 0 such that for 
all 0 < y < y 0 , DfW™ y , < e. Let y 0 > 0 be fixed. For any 0 < y < yo, n e iv™ y and 
corresponding n, we may decompose n as in equation ( 21 ) yielding: 

n i = 2NU*n i + U*w u 1 <i<m (23) 

Since Hndl < a/TA and Until > 1 - yo, it follows that llwdl < 2N^]2yo. Applying con¬ 
straints Ci through C 4 to (23) results in 


Ci: 

llndl 2 + ( nu\ k m ) < 1 + Odlwdl); 

1 < i < m, 


1 < k <N 


llwdl 2 - {nu A k Ui) < 1 + O(HiVill) 

1 < i < m, 


1 < k< N 

C 2 : 

|(n i ,A fe n j ) + (nu\- k nj)\ < CKIM) 

l < i < j < 

m, 

1 < k<N 

C 3 : 

|(ni,A fe nj) - (ni,\- k nj)\ < 0(||Wi||) 

1<i<j < 

m, 

1 < k < N 

C 4 : 

\(n u nj)\ < O(Hwill) 

1<i<j < 

m 



We may thus conclude that there exists a constant C, independent of y 0 , such that n e 
"Wn.cyyo- Since n e I 2 , there exists a T eZ + such that ||rti - II ^ 2yo. Define 


n := (ni. n m ) 

e , Prnt 


nt := (1- 


6a/w ||Prn 


(24) 


It can be shown that for yo < ^, n e ~Wjv, 2 cyjv Furthermore, ||n - n|| < § + 2-Jrhyo. Let 


Cl " 1 “ ( X “ 6Vm) 

R(t) = -C{ l Rn(r), 1 < |t| < JV 

R( 0 ) = Cf 1 (J-J?n(0)) 


(25) 


Then it can be verified that P e _- Cfl . By Lemma 2, for y 0 < 8C ^ : ^ +lj4 , there 

exists signal x el^ such that P x (t) = P(t) for t e [-JV, AT]. 

Define 


Then 


d = n + CiA jV+r+1 x 


(26) 


RAt) = Rn(r) + CiR x (t) (27) 

= 0 t * 0 

= I T = 0 
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Thus d g Ty" 1 . Furthermore, 


||n - d\ 

|| < ||n - n|| + ||h - d|| 


< | + 2^/myo + | 
o 8 

Choosing 


y 0 = min j 

1 Cl e 

1C' 8 C 2 m 2 (N+ l) 4 ’ 6Vm 

gives the required results. 



( 28 ) 


(29) 

□ 


5.3 Image Representations and Implicit Constraints 

We are now in a position to provide a solution to (10) by combining the results of Sections 
5.1 and 5.2. The following theorem outlines how we may replace set 'WJf with TV(f. The 
idea is to penalize disturbance n € TV™ in such a way that the worst case error will occur 
when n is as large as possible; this in effect will force n to be in TV™: 


Theorem 2 Given Mi and M n e Hdfthen 

sup \\Mil + M n n\\ 2 < 1 
le Bl 2 ,ne"Vr™ 

if and only if there exists k 0 such that Vfc > fco. 


sup 

ZeBl 2 ,ne1y™ 


kn 

Mil + M n n 


< mk 2 + 1 


(30) 


(31) 


Proof: 

(31) =j» (30): Since TV™ c It™, it is clear that 


sup sup 

ZeBl 2 nstvy 


k 0 n 

Mil + M n n 


< mkl + 1 


Furthermore, Vn e TV™, |[fco?i|| 2 = mfeg. Thus (30) is implied, as required. 

(30) => (31): Assume (30) is satisfied. Then, by continuity of at y = 0, 

3yo > 0 such that 

sup sup \\Mil + M n n\\ 2 < Co < 1 

leBl 2 

for some constant Co. Furthermore, 

sup sup IIMil+ M n n|| 2 < (IlMdloo + v'?n||M n || 00 ) 2 :=Ci 

ZeBl 2 
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Suppose (31) is not satisfied. Then, Vfc, 3n e It™ and l e BI 2 such that 


1+Cq 

||n|| 2 >m + — 


\\Mil + M n n\\ 2 
k 2 


In particular, let k = Then it follows that n e 4V™ yo . Furthermore, since \\n 

\\Mil + M n n\\ 2 > > Co 

which contradicts our assumption, as required. 


< m, 


□ 


The following corollary may be applied to controller synthesis: 
Corollary 1 Let Go be given. For any controller K, let 

M{K) = [ Mi(K) M n (K ) ] :=G 0 * K 
Then there exists a K such that 

sup \\Mi(K)l + M n (K)n\\ 2 <1 

/eBl 2 ,ne'W'“ 


if and only if there exist K and ko such that Vfc > ko, 


sup 

ieBl 2 ,ne^V 


_1_ 

1 + mk 2 


kVn 

MdK)l + M n (K)Vn 


< 1 


(32) 

(33) 


(34) 


The proof of the above follows from Theorem 2 and lY™ = {n\n = Vn, n e JF}. Next, 
note that constraints Ci to C 4 , which characterize FT, are equivalent to the following 
constraints: 


Cj : lint,kl| 2 < 1, llnull 2 < 1 

C 2 : IInu + n jt k II 2 + llnull 2 + ||%, fe || 2 = \\fn,k + %ull 2 + llnull 2 + llnull 2 
C 3 : l|n fik + % fc || 2 + ||null 2 + llnull 2 = ||hj, fc + null 2 + IInull 2 + ll%,fcll 2 
C 4 : IInu + null 2 + 11%,* + II 2 = llnull 2 + llnull 2 !! + 11%,fell 2 + 11%,fell 2 


1 < i < m, 1 < k < N 
1 <i<j<m, 1 < k < N 

1 < i < j < m, 1 < k < N 
1 < i < j < m, k = 1 


Furthermore, l e BI 2 is equivalent to ||i|| 2 < 1. Thus n <e N and l e BI 2 can be captured 
in a form consistent with 3-f in (11). Finally, let a = , 1 2 . Then by setting 


G = 


0 kaV 0 
aG l n aG^V <xG l2 
G\ 1 G n n V G 22 


(34) may be solved by the method presented in Section 5.1. 


( 35 ) 
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6 Computation 


A controller K may be found which verifies (34) using the methods of Section 5.1 for 
a given k ; thus one would choose k, and synthesize a controller. In order to solve the 
original problem of (33) and (10), however, we need to ensure that k > ko. It is possible 
to find a lower bound for ko given the open loop system G, N, and how closely we want to 
approximate the optimal solution. This bound, however, will more likely than not be too 
conservative to be of any use. In practice, k should be chosen as large as the numerical 
algorithm allows. This is a topic for future research. 

For a given N and m, it is also worth noting how much constraints Ci through C 4 
cost, in terms of the number of constraints (which is linearly related to the number of 
decision variables required): 



COST 

Cl 

2 mN 

C 2 

\m{m - 1 )N 

c 3 

\m(m - 1 )N 

C 4 

\m(m - 1 ) 

TOTAL 

m(2N+ (m- 1)(N + §)) 


Table 1: Cost of Constraints 

Thus the growth is linear in N and quadratic in m. In order to keep the computational 
complexity down, constraints C 2 through C 4 may be omitted, with the result being that 
each component of n will tend to be a white noise signal, but may be correlated to other 
components. The particular nature of the problem will dictate how conservative this 
omission will be. 

In addition, G consists of an augmented version of Go. In particular, V results in an 
extra mJV number of states. Since the number of decision variables grows as the square 
of the number of states of the plant, quadratic growth is unavoidable, both in N and in 

m. 


7 Application to Robust 5/2 Synthesis: Robust Disturbance Rejec¬ 
tion 

We will show the types of robustness problems that may be solved using the machinery 
developed in this paper. Consider the setup of Figure 2. Given P, it is required to design 
K such that disturbances rid, along with measurement errors n e , have a small effect 
on plant output e. The plant is subject to multiplicative, unstructured uncertainty A, 
with associated weight W t . Since na and n e are disturbances, one may know more about 
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Figure 2: Robust Disturbance Rejection 


their spectral content than that they are filtered, but otherwise arbitrary, h signals. In 
many cases, a good model for disturbances is colored white noise, or filtered white noise. 
One may therefore want to model n d and n e as white noise signals, filtered by W )lA and 
W ne , respectively. An example where this may arise is n e being sensor noise and rid an 
impulsive (or known) disturbance (it is also straightforward to consider the case where 
rid is an arbitrary disturbance in BI 2 , etc.). 

The relevant equation describing the system is 

e = SPW t Ae+SPWn d n d + TWn e n e (36) 

= (I- SPW t Ar l (SPW nd n d + TW ne n e ) 

:= Min d + M2n e 


where 


5 := ( I + PK r 1 (37) 

T := PKd + PK)- 1 

The Robust Disturbance Rejection problem is as follows: 

Robust Disturbance Rejection 

For a given N, find an internally stabilizing K such that 

sup sup \\Min d + M 2 n e \\ < 1 (38) 

AeBA (n rf , nJeAVjJ 1 

The above may be converted to a condition which does not involve A: 

Theorem 3 K solves the Robust Disturbance Rejection problem iffK is an internally sta¬ 
bilizing controller and 

sup sup || SPW t l + SPW nA n d + TW ne n e \\ < 1 (39) 

ie Bl 2 (n d ,n e )<=lVff 
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The proof of the above is equivalent to the one in [3], where and n e are assumed to 
be arbitrary l 2 disturbances (the idea is that since A is an arbitrary contractive map, Ae 
may be replaced by l). The above can now be cast into the framework of Section 4. In 
general, many robust synthesis problems may be solved using this technique (see [3]); 
the only restriction is that the uncertainty (possibly more than one block) have e as their 
input. 


8 Example 


We present a simple example for which the solution is known. In particular, we will use 
the machinery developed thus far to synthesize an optimal J-C 2 controller for a given 
plant. We are not proposing this method for synthesizing optimal 3-C 2 controllers, since 
exact solutions exist, but rather wish to explore the properties of our algorithm on a 
simple example. 

Consider the following unstable, non-minimum phase, SISO plant: 


P = 



(40) 


The goal is to minimize the Jf 2 norm of the sensitivity function, 5 = (1 + PK) ' . The 
generalized plant Go for this problem is 


Go = 


1 

1 


-P 

-P 


(41) 


The following table summarizes the results obtained using standard synthesis methods 
(see [4]): 



3-Coc Analysis 

!Pf 2 Analysis 

Optimal Controller 

3-foo Synthesis 

1.1429 

1.1429 

v - 5 

A - 12 

J-L 2 Synthesis 

1.7143 

0.9897 

K = 1 + 0.3889A 


Table 2: Exact synthesis results 

We then repeated the synthesis for k = 0,1,3,10,100 with N fixed at 4, and for 
N = 0,1,2,3,4 with k fixed at 100. The resulting closed loop responses are depicted in 
Figure 3, along with the optimal 3P 2 and 3-foo results for comparison purposes. 

For the first plot, the flat response is that of the optimal 3~C X closed loop, while the 
response with the highest peak is that of the optimal J-f 2 closed loop. The five values 
in between correspond to different values of k, ascending values of k corresponding to 
ascending values of the 3-C x norm. 
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Closed Loop Responses, 


Closed Loop Responses, k=1QO 




Figure 3: Synthesis with constraints 

In the second plot, the 3f x and 3~C 2 designs are also included for comparison pur¬ 
poses. Ascending values of N correspond to ascending values of the 3f 00 norm. Note 
that only six curves seem to be present in the second plot, since the N = 0, k - 100 
design is virtually identical to the 3{ 00 design, as expected. 

The N = 4, k = 100 controller was 

- 1.0079 + 0.3953A 

*' 1 + 0.0306A (42) 

which is very close to the optimal J-C 2 controller. Furthermore, the closed loop 3-f 2 norm 
for this design was 0.9899, and the 3~C X norm was 1.690, again extremely close to the 
optimal 3-C 2 design. Note, however, that the N = 4 design isn’t necessarily the best one. 
For example, the N = 1 design has a closed loop 3d 2 norm of 1.021 and 3d «, norm of 
1.377; As a percentage, the N = 1 design has 3-1 2 norm which is approximately 3% larger 
than the N = 4 design, but 3f x norm which is 23% smaller. The point is that most of the 
reduction in the 3l 2 norm occurred with only one constraint; pushing harder to reduce 
the 3 -C 2 norm only serves to increase the 3foo norm. 

There is an interesting connection between the tools developed in this paper when 
applied to standard 3{ 2 optimization and known results on 3f 2 synthesis. In particular, it 
is known that synthesizing the optimal 3C 2 controller for a given system is equivalent to 
synthesizing the optimal 3f< x controller for the same system multiplied by an appropriate 
weight. Thus our method can be thought of as searching for a feasible controller and 
the correct weighting function simultaneously. This interpretation is also valid for the 
general case, ie., when there is also a signal l. 
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9 Conclusions 

The tools presented in this paper allow for controller synthesis when the disturbances are 
a combination of arbitrary, norm bounded l 2 signals, and signals in set Tk™. For a fixed 
N, the solution takes the form of an LMI, for which good numerical packages exist. To 
solve the problem where noise signal n is perfectly white (ie., satisfies an arbitrarily large 
number of correlation constraints), N will have to be arbitrarily large; this seems to be 
unavoidable, however, since the optimal controller will probably be i nf inite dimensional. 
The strength of this approach, however, is that designing for a fixed N gives us both 
sufficient and necessary conditions for the existence of a controller. By inspecting the 
nature of 'W™ and the resulting closed loop system, the design engineer may then decide 
whether the solution is acceptable, or a larger value of N is required; this is the case in 
the example presented, where it isn’t clear whether the N - 4 design is better than the 
N = 1 design. 

This relates to the view that optimization techniques should only be used as design 
tools, and should be flexible enough to allow the design engineer to customize the results 
to the particular problem at hand; for real problems, the objective is almost never to 
minimize the closed loop or J-f 2 norms (or anything in between, for that matter). 
The best one can do is to provide flexible enough tools to automate some of the design 
process, and let the engineer do the rest. This is the case with standard 3-Coo theory, where 
the problem of optimal control design is reduced to that of finding appropriate weights; 
the design becomes iterative, where weight selection is determined by the previous closed 
loop system. The tools provided in this paper can thus be thought of as a method to 
automate some of the weight selection process, and as a result, allow for more complex 
designs to be performed. 
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Appendix 

Proof of Lemma 2 (see [16]): Let 



' 35?T*(0) 

&*<!> 

R(N) 


Q = 

£*(-1) 

5:1 

/“N 

O 

\R{N- 1) 

g |m(iV+l)xm(N+l) 


! R(-N) \R(-N+ 1) 
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By Gergorshin’s Circle Theorem (see [20]), we may bound each eigenvalue of Q from 
below by 


1 - m y 
N + 


N 


t- - Z m y > 

1 k= 1 


1 - ym(N + l) 2 
N + l 


(44) 


Since y < ^(d+i) 2 > Q > °- Thus there exists P > 0 such that P 2 = Q. Let 


P:-- 


P 

0 


|^2m(JV+l)xm(JV+l) 


(45) 


i.e., we have added m{N + 1) rows of zeros to P. Now define x(t) as follows 

x(0) 


x (2m{N + l) 2 ) 


:= vec ( P ) 


(46) 


i.e., formed by stacking the columns of P into one long vector. We claim that R x ( t) = 
R(t ) for |r) < JV. Partition P as 


Poo 

Pno 


PON 

Pnn 


(47) 


which implies 


N 


R(T) = (N+l-\T\)XPikPkj, j = i + r 

fc =0 


Further partitioning each Pit as 

Pik = [ (Pik) 1 

and using the fact that Pkj = P] k , we have 


( Pik)m ] 


(48) 


(49) 


N m 


R{t) = (N+ 1 - |T |) X I ](Pik)l(Pjk)l j=i + T 

k=0 1=1 


(50) 
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It thus follows from the definition of x(t) in (46), for 0 < t < AT, 


N m N-t 

Rx(t) = X S 2 (Pij)l(Pi+T,j) 
j=0 1=1 i= 0 
JV-t N m 

= X S X (Pij'lltPi+Tj) 

i-0 j=0l=l 
N-t , 

- Zw 7T^ r(t) 

i=0 

= R{ t) 


( 51 ) 


A similar argument holds for negative t. 


□ 
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